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Abstract 

The Black-Scholes implied volatility skew at the money of SPX options 
is known to obey a power law with respect to the time-to-maturity. We con¬ 
struct a model of the underlying asset price process which is dynamically 
consistent to the power law. The volatility process of the model is driven 
by a fractional Brownian motion with Hurst parameter less than half. The 
fractional Brownian motion is correlated with a Brownian motion which 
drives the asset price process. We derive an asymptotic expansion of the 
implied volatility as the time-to-maturity tends to zero. For this purpose 
we introduce a new approach to validate such an expansion, which enables 
us to treat more general models than in the literature. The local-stochastic 
volatility model is treated as well under an essentially minimal regularity 
condition in order to show such a standard model cannot be dynamically 
consistent to the power law. 


1 Introduction 

The Black-Scholes implied volatility is a nonlinearly transformed price of a 
call or put option in such a way that the transformed value does not depend 
on the strike price and the maturity of the option only if the underlying asset 
price is log-normally distributed under the pricing measure. As a function of 
strike price and maturity, the implied volatilities form a surface which visually 
characterizes the marginal distributions of the underlying asset price under the 
pricing measure. In particular, it gives an idea how the price dynamics deviates 
from the Black-Scholes model. Its overall level tells how the underlying asset is 
risky. The implied volatility surface of SPX option prices is usually not flat and 
typically exhibits downward slope and convexity as a function of the log-strike 
price. Further, it has been reported that the slope around at-the-money obeys 
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a power law with respect to the time-to-maturity: see Alos et al. 0, Fouque et 
al. G3, Gatheral et al. Ifl2l . Denoting by at(k, 6) the implied volatility at time t 
with log-moneyness k and time-to-maturity 6, the power law can be formulated 
as 


a t (y/ez,e)-ot(yfe c, e) 
V6(z - Q 


A t d H ~ 1/2 a.s. 


(1) 


as 6 —> 0 for z t ( and t > 0, where H e (0,1 /2) and A is a stochastic process. 
The left hand side is essentially dkOt(0,6). The above finite difference form is 
more relevant here because only a finite number of strike prices are listed in 
real markets. The aim of this study is to construct a model which yields {1}. 

Needless to say, a financial practice needs a model being consistent to the 
implied volatility surface. A popular approach for practitioners is to model the 
underlying asset price S under the pricing measure by a local volatility model 


dS f = S t v(S t , t)dB t/ 


( 2 ) 


where v is a Borel function and B is a standard Brownian motion. As shown 
by Dupire |5|, for any arbitrage-free set of vanilla option prices, there exists a 
function v such that each of the given option prices coincides with the theoretical 
no-arbitrage price under (0 for the corresponding payoff. The procedure of 
finding such a function v given a set of market prices, called the calibration, 
has been the first step to price exotic derivatives without providing a static 
arbitrage opportunity. 

A more important practice is to hedge an option portfolio. The local volatil¬ 
ity model is not satisfactory for this purpose due to the lack of dynamic consis¬ 
tency; the calibration to market prices at different times usually gives different 
functions as v. As a result, an hedging strategy under a model calibrated at 
time f is outdated at time s > t. This simply implies that the underlying asset 
price S in fact does not satisfy (0. The hedging error is then out of control, at 
least from theoretical point of view. 

The necessity of the re-calibration can be deduced from the lack of dynamic 
consistency to the power law. As shown in Section 2, under a local-stochastic 
volatility model extending 0 : 


dS f = S t v(S t ,Y t ,t)dB t , 

k 

d Y\ = b l (S tr Yf, t)dt + ^ c)(S t , Y„ t)dwj, i = 
i =i 


where W = (W 1 ,..., W ,: ) is a k-dimensional standard Brownian motion with 
d(B, W) t = p‘(St, Yt, f)df, we have 


a t (Vez,e)-ot(Vec r e) 
V6(Z - Q 


-y 

- {Sfd s r?(S f , Y f , t) + ((cp) ■ V y log v)(S t , Y t , f)} a.s. 

(3) 
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as 9 —* 0 for all t and z + C This result partially extends Medvedev and Scail- 
letQi and Osajima (20) . The point here is that © holds under an essentially 
minimal regularity condition. This implies that the local volatility model |2]l 
needs a volatility function v(s, t) which is singular at (s, f) = (St, t) in order to 
be consistent to the power law at time t. To be dynamically consistent, the 
model needs a volatility function v(s, t) which is singular everywhere and this 
is nonsense. 

As an application of general theories, Alos et al. f2| and Fukasawa (lOl 
treated volatility processes driven by a fractional Brownian motion to find that 
the term structure CD> at fixed time t follows under those specific stochastic 
volatility models. The Hurst parameter H of the fractional Brownian motion 
has to be chosen from (0,1/2) to match Therefore the volatility is not a 
process of long memory. See Gatheral et al. |fI2| for an empirical work which 
suggests that the volatility appears in fact a fractional Brownian motion with 
H e (0,1/2). The models of Alos et al. J2] and Fukasawa flOl are however not 
dynamically consistent to (HI in the sense that the models have to depend on 
t to yield (JTj) . Therefore it suffers from the same drawback as local volatility 
models do. Due to the fact that the fractional Brownian motion is not Markov, 
the construction of a dynamically consistent model is not a trivial exercise. In 
Section 3, we use a representation of a fractional Brownian motion given by 
Muraviev H8l to solve the problem and show that the constructed model in 
fact yields CL} for all f. More precisely, we have ((TJ) for all t with 

At = cdlogc>(Y f ) (4) 

under 

d S t = S t v(Yt)dB t , 

r t (5) 

Yt = Y 0 + b(Y u )du + W?, 

Jo 

where W H is a fractional Brownian motion with Hurst parameter H e (0,1 /2) 
with correlation 

E[(B t+ g - B t )(W? +g - Wj 1 )\Tt\ = c' 6 h+1/2 (6) 

and c, c' are constants. The potential usefulness of the Muraviev representation 
in finance was discussed by Novikov |[T9| . 

Now we explain how our asymptotic analysis is related to other approaches 
in the literature. There are two major categories. The first one is based on a 
perturbation from the Black-Scholes model. The idea is to introduce an artifi¬ 
cial perturbation parameter in, say, a stochastic volatility model in such a way 
that the model converges to the Black-Scholes model in a suitable sense as the 
perturbation parameter tends to zero. The implied volatility then converges 
to the volatility parameter of the limit Black-Scholes model. An asymptotic 
expansion of the implied volatility is then derived around the limit volatility 
parameter. The small vol-of-vol expansion by Lewis JT5], the singular pertur¬ 
bation (fast mean reverting) expansion by Fouque et al [7] and the multi-scale 
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expansion by Fouque et al j8j belong to this category and can be verified in a 
unified manner as in Fukasawa [10j|. The second category is based on the short- 
time behavior of the underlying asset price process. There are again two major 
approaches within this category. The first one considers the implied volatility 
in the original scale of the strike price. The large deviation principle and the 
heat kernel expansion are explicitly or implicitly underlying this approach. See 
e.g., Berestycki et al 0], Osajima [20] 12T1 , Henry-Labordere fl3l . Pham [22| , 
Forde and Jacquier [6], Gatheral et al (Til and Armstrong et al. J3). The sec¬ 
ond one rescales the strike price to get a high resolution around at-the-money. 
See Yoshida [23], Kunitomo and Takahashi Ifl4l , Medvedev and Scaillet 06] , 
Osajima I20l and Mijatovic and Tankov [17]. This last approach is the most 
relevant here because we are considering the term structure of the implied 
volatility around at-the-money. As already mentioned, I0 partially extends 
Medvedev and Scaillet [16] and Osajima [201 . The former is based on a formal 
expansion of the PDE that the implied volatility satisfies. The latter, as well as 
Yoshida 1251 . Kunitomo and Takahashi [14] , is based on the Watanabe theory 
of the Malliavin calculus. Our new method requires less regularity conditions 
and is effective to derive (JT} under 0. We do not take the jumps of the asset 
price process into account because Mijatovic and Tankov [17] has already con¬ 
sidered an exponential Levy model to find that the implied volatility behaves 
differently from ((TJ. 

We conclude this section with an additional remark on the jumps. By a 
seminal work by Ait-Sahalia and Jacod [|T|, it has been widely believed that 
continuous asset price models are rejected by statistical testing. In fact, Ait- 
Sahalia and Jacod [Tj and other related studies have always assumed that the 
volatility is an Ito process. In other words, what has been rejected is only a 
continuous asset price model with volatility being an Ito process. The rough 
fractional volatility in (O is not an Ito process. It remains for future research to 
develop a statistical theory for such a model based on high frequency data. In 
this direction, as already mentioned, Gatheral et al. [T2I found that modeling 
the volatility with W H , H e (0,1/2) is consistent to a scaling law observed in 
high frequency volatility times series. 


2 The local stochastic volatility model 

Here we study the short-term behavior of at-the-money skew under a regular 
local stochastic volatility model extending 0. Let (Q, T, P, {Tt\t>o) be an fil¬ 
tered probability space satisfying the usual conditions. We suppose a Markov 
structure under the pricing measure: 


dSf = Stv(S t , Y t , f)dB ( , 

k 

d Y\ = b\S t , Y t , t)dt + ^ c)(St, Yt, t)dW j t , i = 1. d, 

i =i 


4 


where B is an {^[‘Standard Brownian motion and W = (W 1 ,..., is a k- 

dimensional {^[-standard Brownian motion with d(B, W')t = p'(St,Yt,t)dt. The 
continuous functions v, b‘, c l . and p 1 are defined on (0, oo) x X [0, oo). Denote 


1=1 
n d\ n — 


1=1 


•£<& t-t.#' v »=(^-.^ 


d_ 

dy a 


and i] = (tj, ..., if), a = [rf ! |' j=| . Here we work under the following regularity 
conditions: 


1. v(s, y, t) is positive, bounded in s and of linear growth in y 

2. b\s, y, t) and c' (s, y, t) are of linear growth in (s, y). 

3. v(s, y, t) is continuously differentiable in (s, y) and that there exists k e N 
such that 

\d s v(s,y,t)\ + \V v v{s,y,t)\ 
sup -—-< oo. 

s>0,yeR ,i ,f>0 1 T s 


4. v(s, y, t) is locally H-Holder continuous in t with H > 1/2. 


Theorem 1 For any zeR and t > 0, 


as 6 —» 0, where 

e Vez _ i 

A =- v t = v(S t , Y t , t) 

Vd 

and 

at = ^(v(S t ,Y t ,t) + S t d s v(St,Y t ,t) + (i) ■'V y logv)(St,Y t ,t)). 


Proof: Since (S, Y) is a (time-inhomogeneous) Markov process, we can and do 
assume t = 0 and To is trivial without loss of generality. Define the rescaled 
processes X e and Y e by X® = 0 _1/2 (Se„ - So)/So and Y® = 9~ 1/2 (Yg u - Yo). The 
rest consists of two steps. 


Step 1) Here we show that (X®, Y®) is uniformly integrable in 6 and converges 
in law to, say, (X®, Yf) which is normally distributed with 

E[X°] = 0, 

E[Y°f] = 0, 

E[\X° u \ 2 ] = v(S 0 ,Y 0 ,0)\ 

E[X°X' i ] = v(S 0l Y 0r O) n \S 0l Y 0l O)u, 

E[Y°,fY° u ’] = a'i(S 0 ,Y 0 ,0)u 
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as 9 —» 0 for each u > 0. Under the regularity conditions, an application of 
Gronwall's lemma gives that 

E[sup |S„| 2 ] + £[sup |Y„| 2 ] < oo 

0 <u<t 0 <u<t 

for each t > 0. It follows then that 

-2 


and so. 


£[ix; 


lim E[SgV(Sg, Yg, 9) \ = Sfr(So, Y o ,0) z 
6—>0 


»I 2 1 = V f E l s lr v ( s 8r , Yer, 9r) 2 ]dr v(S 0 , Y 0 , 0) 2 » 

Jo 


as 6 —» 0 for each u > 0. In particular, we have that X® is uniformly integrable 
in 0 for each u. Let M® be the local martingale part of Y®. Since 


e[|m; 


S|2i 


d d 

= Y mHSeu / Y e „, 9u)\ 2 ]du -» Y l«"( s o> y o> »)l 2 » 

;=i i =i 


as 6 —» 0, the vector M® is uniformly integrable in 6 for each u. Further, the 
absolutely continuous part of Y® converges to 0 in E 2 . Note also that 

r*U 

( X®>„ = Sq 2 S 2 gr v(S d r, Yer, 9r) 2 dr -» p(S 0 , Y 0 ,0) 2 u, 

Jo 

r»U 

<X®, Y®)„ = Sq 1 SeMSer, Yg r , 9r)T } (S er , Yer, 9r)dr -» z>(S 0 , Y 0 ,0)/](S 0 , Y 0 , 0 )h, 
Jo 

l 

a(S gr , Ye r , 9r)dr -> a(S 0 , Y 0 ,0)w 


<Y®, Y®) 


r u 

Jo 


as 9 —> 0. Therefore the convergence in law follows from the martingale central 
limit theorem. 


Step 2) Define p by 


p(x,») = (A - x)® 


A - x 


U>o Vl - » 

This solves the partial differential equation 


+ fo 


Vi - nip 


A - x 


fo 


vr 


Note also that 


d x p(x, u) = -O 
<? 2 p(x, u ) = 


j = °/ 
p(x,l) = (A - x)+. 

A - x 


Vo 


VT“ 


Vo 




U Wo 


A - x 


iVT^ 
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Let f(x, y, t) = (1 + x) 2 u(So(l + x), Yo + y, t) 2 . Then, by Ito's formula, 
E[(S 0 e^ z - Sg), 


s 0 Ve 


=E[( A-X“) + ] 


=p( 0, 0 )+^f g E[d 2 x p(x e u , u){f{ Vex®, Ver®,e«) - ? 2 )]du. 


Since 


d 2 p(x, u) 


-.(p 


Z - X 


Vo V1 - H \ l? 0 Vl - M 
as e —> 0, applying the result from Step 1, we have 

e- 1 / 2 E[e 2 P (x®, »)(/(Vex®, Ver®, e«) - /(o, Ver®, e«))] 

= E[e 2 P (x®, M )x® f d x f(A Vex®, Ver®, e«)dA] 

Jo 


d x f (0,0,0)- 


v 0 


VT 




vr 


Since 


■du < cc, 

Vl-u 

the dominated convergence theorem gives that 


u L \»o 


r 1 — 

Jo Vl - ; 




^f o E[d 2 x p(x e u , u)(f( Vex®, Vey®, e«) - /(o, Ver®,e M ))]du 
= |(o 0 + S 0 e s p(s 0 , Y 0 , 0)) Ve</> + o( Ve). 


Similarly, 


e- 1 / 2 E[d 2 P (x®, H )(/(o, VeY®, eu) - /(o, o, e«))] 

= £[d 2 p(X,®, U ) f Y®-V y /(0,AVeY®,e»)dA] 
Jo 


Vo 


vrv 


r z — X° 

<H A=== 1 Y» ' V y /(0,0,0) 


M L Wo 


,vr 


and so. 


= (^) T7(S 0 , Yo, 0) • V y /(0,0,0), 

\ f o E[d 2 x p(x° u , u)(f( o, VeY®, du) - /( o, o, e«))]d« 
= 2 ^ (f? ' V ^ )(S °' Y °' 0) ^ (^) + 0( ^ 


(7) 
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Finally, since /(0,0, •) is locally H-Holder continuous with H > 1/2, 


Liri ii 

which completes the proof. 



E[cfy( x i w)(/( 0, 0 , Oil) - /(0, 0 , 0))]d« = o( Vfj), 


//// 


The Black-Scholes price at time t of a put option with time-to-maturity 0 and 
strike price K = S t e k is by definition 


P t (k,e,a) = S/®(-rf 2 ) - S f <D(-di), 


where ® is the standard normal distribution function. 



oVd 


This is an increasing function of the parameter a, which is called the volatility. 
The Black-Scholes implied volatility Ot(k, 0) at time t is defined through the 
equation 

P t (k,6,o t (k,6)) = P t (k,6), (8) 

where Pt(k, 0) is the corresponding price of the put option. 

Theorem 2 For any zeR and t > 0, 


a t (V0z,0) = v(S f , Y t/ 1) 



as 0 —> 0. 


Proof: In light of Theorem [lj it suffices to show that 



implies 



Using that A = z + V0z 2 /2 + o( \fd) and cf)'(x) = -x<p(x), we have 
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On the other hand, by definition, 

Pt( yfdz, 9, a) AA , JX . ®(-d 2 )-®(-d!) 

- — - = AO(-fl 2 ) +- ■= - 

s,Vd Vd 

= AO I - j + cp ^ Vdz + o(p j + 0 ( Vd). 

It is then easy to see Ot( V0z, 8) = Vt + o(l). Put /3(d) = a t { Vdz, 8) - v t . Then, 

P f ( Vdz, d, v t + yS(d)) 

S f Vd 

= AO (|) + 4>(f)|Vdz + (o f + /3(0))</> (|) + o( Vd) + o(/8(0)) 

and so, 

at Vd = y Vdz + /3(d) + o( Vd) + o(p(8)). 

Therefore, 

d“ 1/2 /3(d) -> a t - ~z, 

which completes the proof. //// 

We have seen that if z ^ 

ff t (Vdz,d)-CTf(C,d) 1 J c j v ,s , , „ , wc v 

-p-» - S t d s v(S t , Yt, t) + (i] ■ V y log v)(S f , Y f , t) a.s. 

Vd( z-Q 

as d —» 0 under a regular local-stochastic volatility model. This means that the 
regular models cannot explain the empirically observed term structure of the 
implied volatility skew (T]i. As a result, a forcible calibration at time t leads to 
a function v which has a kind of singularity at (s, y, t) = (St, Y tl t). Therefore, 
it has to be singular everywhere in order to have (Q} for all t. It is hopeless 
and even nonsense to have such a singular function by any practical method of 
calibration. 


3 Rough stochastic volatility model 

Here, we consider a model with stochastic volatility being rough. More pre¬ 
cisely, we model the volatility to be driven by a fractional Brownian motion 
with Hurst parameter H e (0,1/2). We show that the term structure ijl) follows 
when the fractional Brownian motion is correlated with the Brownian motion 
that drives the underlying asset price process. Let (LI, T, P, \'Ft\t>o) be an filtered 
probability space satisfying the usual conditions. The underlying asset price 
model under the pricing measure is 

d S t = Stv(St,Yt,t)dB t , 

Y t = Y 0 + f b(Y u )du + W f H , 

Jo 
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where B is an {^[-standard Brownian motion, W H is a fractional Brownian 
motion with Hurst parameter H e (0,1/2). Here we work under the following 
regularity conditions: 

1. v(s, y, t) is positive, bounded in s and of linear growth in y. 

2. v(s, y, t) is continuously differentiable in (s, y) and that there exists k e INf 
such that 

\d s v(s,y,t)\ + \d y v(s,y,t)\ 
sup ---< oo. 

s>0,y€R,f>0 1 + S 


3. for each (s, y), there exists e > 0 such that v(s, y, t) is locally (H + e)-Holder 
continuous in t. 


4. for each s, there exists e > 0 such that d y v(s, y, t) is locally e-Holder con¬ 
tinuous in (y, t). 

5. b is Lipschitz continuous. 

To introduce the correlation between B and W H , we adopt a representation of 
the fractional Brownian motion given by Muraviev ||T8| : 

w? = CH f o ^(Zf-zjw, 

zf = f e“^ (t “ s) dW s , 

oo 

where Ch > 0 is a constant and W is an {^(-standard Brownian motion. We 
assume d(B, W) f = p(Y t )dt with a continuous function p with \p\ < 1. 


Theorem 3 For any z e 1R and t > 0, 

-S„ 9 ),r,l =to (A| afgH M' + as 

.<?, \v t j \Vt) 


as 0 —> 0, where 

e Vez _ i 

A =-—, Vt = v(S t , Y f , t) 

Ve 

and 


a t = (i/2 fi +HX3/2?H) 2p(Yf)<9y l0§ U(Sf ' Yt ' + F t d y v ( St ’ Yi ’ F >' 

/~*oo 

F e = I ^-3/2-H (1 _ p _ e-^o-^z^dp. 

Jo 

Remark 1 F® is afunctional of{ W s }_oo< s <f- T/ze lazv ofF ® does not depend on 9. 
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From Theorem[3l by the same argument as in the proof of Theorem[2j we obtain 
the main result of this paper: 

Theorem 4 For any zeE and t > 0, 

o t ( Vdz, 6) = v(S t , Y it t) + a e t 6 H + o(d H ) a.s. 

as 6 —> 0, where af is defined by ©. In particular, if z + f, 

as 6 —> 0. When p is constant and v(s, y, t) = v(y), we have 0 with 0 and 

We start with some lemmas. 

Lemma 1 For every f,seR, 

/~%00 /-%00 

I |3- 1/2 - H |zf - Z^IdjS < co, I p-V2-H\ e -iKt-s) _ i\\2f s \dp < co a.s. 

Jo Jo 

and 

/-•co 

Wf = Wf + c H I r 1/2 ~ H (Zf - z s )d/3 
Jo 

/-»oo 

= Wf + c H |S“ 1/2 “ H (e“^“ s) - l)Zfd|3 

Jo 

X oo r*t 

p-V2~ H J e -m~u) A w u djS a.s.. 

Proof: Let t > s without loss of generality. By definition, 

Zf = + J e-P (t ~ u) dW u 


and so. 


E[|Zf - Zf | 2 ] < 2 |e^ s - f) - l| 2 E[|zf | 2 ] + 2 J e - 2 W-")d M 
= ^“V (s - f) - 1| 2 + f~\l - e 2 ^ (s_f) ) 

= 2 J 8“ 1 (1 - e _ ^ (f_s) ) 

< 2min{(f - s),/3 -1 }. 


Therefore, 


f 


Loo 

The rest is obvious. 


f |Zf - Zf|d^ 


/-*00 

I jr 1 / 2 - H E[|Zf - zf| 2 ] 1/2 dyS 
Jo 


< oo. 


//// 
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Lemma 2 For all t el R, 0 > 0 and for all Borel function f, 

E[f(S t+ e,Yt + e)\rt] = E[f(S t+ g,Y t+ 0 )\S t ,Y t AZ p t } l}>o ] a.s. 

and 

E[f(S t+e ,Y t+e )\St = s,Y t = y, zf = > 0] 

= E[f(Sg,Yg)\S 0 =s,Y 0 = y,zj = zP,p > 0] a.s.. 

Proof: Since 

Y u =Y t + J b{Y r )dr + Wf - Wf 

=Y t + J b(Y r )dr + c H ^-i /2 -« e -^(»-0zf d^S 

X OO pU 

p- 1/2 ~ H J e -«“- r) dW r dj3 

by LemmaU the conditional law of {Y u } u > t given St, Y t/ {Z’f \p>o is independent 
of Tt- Therefore for any bounded ^-measurable random variable A, 

E[f(S t+ g r Y t+e )A\S t/ Y tr {Z^> 0 ] 

= E[f(S t +g, Y t+e )\S t , Y t , {Z p t }p> 0 ]E[A\S t , Y t , [Z% >0 ] a.s.. 

It follows then that 

E[ f(S t+g , Y, +e )A] = E[E[f(S t+ g, Y t+e )\S t , Y t , {zf } p> 0 ]A], 

which proves the first part. The second part then follows from the stationarity 
of the Brownian increments. //// 

The above lemma shows that all the information on the history of W S/ s < t is 
translated into the spot values z| ! , fi > 0. We may say (i i-» z!’ t is the Laplace 
transform of W S/ s < t. For each fi > 0, Z? is an OU process 

dz\ = -ySZfdt + dW f 

and in particular, a Markov process. Note also that the infinite dimensional 
process 

(S, Y, Zf-, j3 > 0) 

is a time-homogeneous Markov process. 

Proof of Theorem [3j By Lemma |2[ we may and do assume without loss of 
generality t = 0 and the existence of a regular conditional probability measure 
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Po given To- In particular, v = v(So, Yq, 0). Denote by Eg the expectation with 
respect to P 0 . Let X® = d~ 1/2 (Sg u - S 0 )/S 0 , 

r»oo r*Qu 

y® = c H e~ H I jg- 1/2 -» I e -/s(eK-b dW) . djg 

Jo Jo 

and 


YS = 0 


-H 


r*Qll 

Jo 1 


fc(Y r )dr + Y®. 


Then by Lemma|l] we have the decomposition 


J r*0u 

0 


= Yq + ch 


o 


r 

Jo 


You =Y 0 + I b(Y,.)dr + CH I _ i)zgdjS + 0 H Y® 


j8 - l)ZgdjS + 0 H Y®. 


The rest of the proof consists of three steps. 


Step 1) Here we show that (X®, Y®) is uniformly integrable in d under Po and 
converges in law to, say, (X®, Y®) which is normally distributed with 

E[X°] = E[Y° u ] = 0, 

£[|X®| 2 ] = v 2 0 u , 

£[X®Y®] = vop(Y 0 )ch / 2 + u H+1 ' 2 ■ 

as 0 —» 0 for each u > 0. Under the regularity conditions, an application of 
Gronwall's lemma gives that 

E 0 [sup |S„| 2 ] + £ 0 [sup |Y„| 2 ] < oo 

0 <u<t 0 <u<t 


for each t > 0. It follows then that 

r»6u 

d~ H b(Y r )dr —> 0 

Jo 

as 8 —> 0 in L 2 (Po). By the scaling property of W, the law of Y® under Po does 
not depend on 6 . It follows then that (X®, Y®) is uniformly integrable under Po- 
Let W® = 8 ~ 1 / 2 Wg u . Then, we have 


r*oo r*U r*U 

E[W,®Y®|r 0 ] =c H I p~ 1/2 ~ H I e -«»-d dr d J 8 = c H r(l/2-H) (u - r) H " 1/2 dr 
Jo Jo Jo 


=Ch 


m/ 2 -H) H+ll2 

1/2 + H 


13 




Note that B admits a representation 

X U r*U I - 

p(Y r )dW r + ^1 - p(Y,) 2 di< 

with (W, W^) being a 2-dimensional {^[-standard Brownian motion. Then, 

X® - vp(Y 0 )W e u - v - p(Yo) 2 6 ~ 1 / 2 Wg u -» 0 

in L 2 (Pq) as 6 — > 0. It follows that the law of (X®, Y®) under Pg converges weakly 
to the law of (X®, Y,° ( ). 

Step 2) Let g(x, y, z, f) = (1 + x) 2 »(So(l + x), Yo + y + z, t ) 2 and define the function 
p as in the proof of Theorem[l] Then, by Ito's formula, 

£q[(V^ - S 0 ) + ] 

s 0 Ve 

= £o[(A-Xf) + ] 

= p(o, + Eo[<? 2 vP (x®, «)(g( Vex®, e«Y®, hm, eu) - v 2 0 )]du, 


where 


Since 


h(u) = c H 
d 2 x p(x, u) 


r 


P ~V2-H( e -j}u _ l) Z jjd/3 

1 


--<P 


z- X 


Vo V1 - It \ z>0 Vl - w 
as 0 —» 0, applying the result from Step 1, we have 

er 1/2 E 0 [^ P (x®, u)( g ( Vex®, e H Y®, fi(e«), e») - ^o, e H Y®, £(e«), e«))] 

= E 0 [e 2 P (x®, M )x® f d x g(\ Vex®, e H Y®, fc(e«), e«)dA] 

Jo 


d x g( 0 , 0,0,0)- 


VT 


x° 


vr 


x 


= <?xg(0,0,0, 0)—(p (—]. 

v 0 \v 0 J 

Therefore, by 0 and the dominated convergence theorem, 

\ f 1 E 0 [^P(X®, u)(g( Vex®, e H Y®, h( 6 u), eu) - g (o, e H Y®, fi(e M ), e«))]d« 
= o(e H ). 
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Similarly, 


e- H EoK-p(X^, u)(g(0, G H Y°, h(6u), Gu) - g( 0,0, h(6u), 0m))] 
A,, v e ,a P Y®d y g{0,AG H Y e u , h(6u), 6u)dA] 


H 

= E o [0^(X o e , M ) 

Jo 


Vo 


VT^ 


<p 


; - x° 


»0 


vr^ 


^_I y0 


<?yg(0, 0, 0, 0) 


CH- 


r(l/ 2 -H) zm h+1/ 2 /z 
1/2 + H vl 


^ ) p ( s °' v °' 0) a y^(0' 0' 0' o)' 


and so. 


\ f o Eo ld 2 x p(X e u , u)(g(0, V0Y®, h(Gu), Gu) - g( 0 , 0 , h(6u), 6u))]du 

= CH w^mwm zpn),> ’ '° sv{So ' ¥ °' 0>f>H,p (^) + ° (f>H> ' 

Next, observe that 

1 r 1 

- J E 0 [d 2 x p(X°, n){g( 0,0, h(Gu), Gu) - g( 0,0,0,0M))]du 

= ±f E 0 [G 2 p(X°,u)] d z g(0, 0, Ah(Gu), 6u)dAh(6u)du 

_ 1 
” 2 


E 0 [d“.p(X^, m)] -—</> — 


J z g(0, 0, Ah(Gu), Gu)d/\h(Gu)du 

i \mo//jo 

r —</>(—) r 

2 Jo vo Y \v 0 JJo 
+ G H FgGyV(S 0/ Y 0 ,0)(p (—). 

Wo / 

We show in the next step that 

f (e 0 [JJ>(X®,m)] </>(—)) f d z g(0,0,Ah(Gu),Gu)dAh(Gu)du = o(G J 

Jo ' Wo//Jo 


^o v^o // Jo 
(J 2 g(0,0, Ah(Gu), Gu) - d z g{ 0,0,0,O))dA/z(0M)dM 

A' 

Wo/ 


') 

( 10 ) 


and 



(J z g(0,0, Ah(Gu), Gu) - d z g( 0,0,0,O))dA/z(0M)dM = o(0 H ) 


(ID 


in the almost sure sense. Then, the proof is completed by noting that 


1 

2 


f 

Jo 


Eo[dxP(X%,u)(g(0, 


0,0, 0m) - g(0,0,0,0))]dM = o(0 H ) 
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since y(0,0,0, ■) is locally (H + e)-Holder continuous with e > 0. 


Step 3) Here we show ( flOb and (Hil l. Let 

/-»oo 

h(u) = c H f}- ll2 ~ H \e~P u - l||Zg|dj3. 

Jo 

Let us show that 

9~ H+e h(9) -» 0 a.s. (12) 

as 0 —> 0 for any e > 0. Since 


sup \h(6u)\ < h{9) 

0<u<l 

and d z g(0, 0, •) is locally e-Holder continuous, dTTb follows from (El) . As we 
have already seen in the proof of Lemma[lJ E[h(u)] < oo. Further, by changing 
variable as y = 8f3, 


e~ H h(6) = C H 



y -ll2-H\ e -v _ n 



e yv dW 9 


d y. 


which means that the law of 6 H h{9) does not depend on 6. Therefore, 

E[6~ H+£ h(6)] = C6 e 

with a constant C > 0. Let 9 n = n~ 2 ^ e . Then, for any 5 > 0 


£ p(\e~ H+£ h(e n )\ >5 )<cY j ^- 2 <^ 

n =1 n =1 


and so, 

e~ H+£ h(e n ) -> 0 a.s. 

by the Borel-Cantelli lemma. This implies dl2l l because h is a decreasing func¬ 
tion and 6 n /6 n + 1 —» 1. Now, notice that q(x,r) = d 2 p{x,r) solves the partial 
differential equation 

M + \vffin. ^ (0 ' 0) = ^(^) 

and so, by Ito's formula. 


E 0 [dlp(X e u ,u)]-±<pl^) 

Vo V V 0 / 

= \ f E 0 [d$p(X?, r)(g( Vex®, e H Y 9 r , h(9r), dr) - v 2 0 )]dr. 
Repeating the same argument as Step 2, with the aid of ((12) , we have that 


sup 

0<w<l 


; 0 [d 2 x p(X e u ,u)]--$(-) 
v 0 \v 0 / 


= o(9°) 


for any <5 £ (0, H). This and (fl2ll imply d 1 0b . 


//// 
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